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Abstract. We have calculated the emission spectra of hydrogen and sodium atoms in the cool part of 
prominence models which satisfy simultaneously the constraints of radiative transfer, statistical equilibrium 
and charge-particle conservations. 

In the considered range of our model parameters, emission strengths of Hg and Nat D lines increase with 
the temperature and the total number density. Low-pressure models raise the ionization rate highly but yield 
very weak Nal D line intensities, since these model prominences contain small amounts of free electrons 
and sodium atoms which have a deep relation with the formation of sodium lines. We find that sodium 
D lines should be emitted in the high pressure region of prominences, and that their intensities are difficult 
to attain in the cool core of any model prominence with a temperature as low as 4000 K. In order to explain 
consistently the spectral emissions of Ha and Nal D lines observed in quiescent prominences, a total 
number density higher than 4 x 10!! cm ^? and a temperature over 5000 К are required at least in the cool 
part of prominences. 


1. Introduction 


Spectroscopic studies on solar prominences have been carried out theoretically by a 
number of authors (Hirayama, 1964; Kawaguchi, 1964; Poland et al., 1971; Chultem 
and Yakovkin, 1974; Heasley and Mihalas, 1976). Among them, Heasley and Mihalas 
(1976) have achieved a remarkable work on the hydrogen spectrum and the structure 
of quiescent prominences which satisfy self-consistently the constraints of force and 
energy balance as well as of radiative transfer and statistical equilibrium of hydrogen 
atoms. Similar calculations based on the isothermal and isobaric conditions have been 
made by Heasley and Milkey (1976, 1978, 1983) to construct theoretical models of 
quiescent prominences with respect to the temperature, density and helium abundance. 

In contrast with model structure of prominences simplified as one slab with finite 
thickness, Fontenla and Rovira (1983) have treated the prominence as a group of 
numerous slabs characterized by ‘filament’ or ‘thread’ confirmed observationally in the 
earlier works (Dunn, 1960; Engvold, 1978). More recently, Fontenla and Rovira (1985) 
have presented the result of calculation for the hydrogen lines and continuum emitted 
by a quiescent prominence with a number of filamentary structures, using the constraints 
of energy balance and isobaric condition differing little from those of Heasley and 
Mihalas (1976). Going one step forward from the one dimensional transfer computation 
of the atmospheric model with an infinite slab geometry. Heasley (1977) has solved the 
non-LTE line transfer problem in a circularly symmetric, cylindrical region. Two- 
dimensional transfer calculations of resonance lines in prominence models with an 
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uniform temperature (Т, = 8000 К) and given ionization degree (1 < n,/n, < 10) have 
been made by Vial (1982b), with the modified code of Mihalas et al. (1978). 

Interpretation of optically thin line spectra is rather easy without any inclusion of 
radiative transfer processes, even though we are often confronted with the problem 
caused by the lack of atomic data. For the optically thin lines of metals (Na, Mg, Sr *, 
Ва”) observed in prominences, statistical equilibrium solutions have been obtained by 
Landman (1981, 1983, 1984) who used the recent atomic data given by experiments and 
if not available, his own calculations. In a series of three papers (Landman, 1983, 1984, 
1985), he derived the physical conditions in the cool part of prominences such as a low 
ionization rate of n(H * )/n(H 1) ғ 0.07, high gas pressure of P ғ 3-6 dyn сіп ” 2, tem- 
perature of Т ~ 7000 К and etc. Lately, his results on the ionization degree of hydrogen 
and the gas pressure were revised to 0.16 and 1.4 dyn сш” ?, respectively, by the analysis 
of the effects of model Na atom level truncation (Landman 1986). 

Intensive studies on sodium Nai D lines in solar prominences have been made by 
Nikaido and Kawaguchi (1983) who performed a non-LTE transfer calculation based 
on an isothermal slab model with an assumed ionization rate of hydrogen (0.1). They 
found that the intensity ratio of D1to D2 emission lines deviates significantly from their 
theoretical value in some prominences. For this type of prominences, they have shown 
that the sodium line emitting region should be in high gas pressure and at very low 
temperature (4000 K), in order to explain the large optical thickness of observed lines. 

However, their study should be improved in the following areas. First, the assumed 
ionization rate of hydrogen atoms may be too high for the cool part of prominences with 
a temperature as low as 4000 K, even if a thread-like structure of prominences is 
considered for an easy penetration of Lyman continuum radiation. Second, the distribu- 
tion of free electrons inside the slab should be taken into consideration, since the number 
of electrons varies rapidly from the surface hit by EUV radiations to the cool center of 
prominence. Actually, the distribution of electron density exerts a fair influence on the 
formation of NaI D lines via the recombination processes. Third, according to observa- 
tions, comparatively intense emission of NaI D lines is usually accompanied by strong 
hydrogen Balmer lines. Then, prominence models should be constructed to match 
successfully such an observational fact. 

The purpose of this paper is to investigate the physical conditions of solar 
prominences which can explain consistently the observed emission spectra of Ha, Hf, 
and Мат D lines. In order to take a full account of the above mentioned problem, we 
have solved the non-LTE problems of hydrogen and sodium atoms simultaneously. We 
will present the result of our calculation on hydrogen and sodium line emissions for 
prominence models with isothermal and isobaric slab structure. 


2. Method 


2.1. BASIC EQUATIONS 


The emitting region of sodium lines in solar prominences has been examined by using 
a stationary model in which ionization and thermal equilibrium is supposed to be 
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satisfied. We consider for simplicity that the prominence embedded into the corona is 
made up of a cool core with uniform density and temperature inside. The geometry of 
the cool core we consider is an infinite slab with a finite thickness, symmetrical with 
respect to the central plane and illuminated on both sides by the incident radiation 
coming from the photosphere, the chromosphere, the transition region and the corona. 
We will treat the cool core as a target model of prominence to investigate from now in 
our work. 

The basic equations used in our model are as follows: 

Transfer equation: 


02(7,7,)/дт2 =J,- Sy. (1) 
Rate equation: 
л, (X )/dt = F,(J,, п) - DJ, n). (2) 


Charge conservation: 

n, = Xn(X'). (3) 
Particle conservation: 

N =n, + ZN(X). (4) 


Here f,, J,, and 5, are the Eddington factor, mean intensity and source function at 
frequency v, respectively. F, and D, are the total formation and destruction rates of 
atomic species X in the excited state of energy level i, which are generally expressed by 
both radiative and collisional transition processes. n;(X) is an occupation number for 
level i(i = j). The net rate of n;, dn;/dt is, of course, set to zero, by the assumption of 
stationary state. And n,, n(X * ), N(X), and N stand for the electron number, ion number, 
number of species X and total particle number density of the model, respectively. In the 
equation of charge conservation, we ignore the number of free electrons which come 
from the second ionization stage of the species. 

After the transfer equation is transformed to difference equations, all equations are 
linearized for the perturbations of unknown quantities. And then, we solve simul- 
taneously the full set of above equations, in a similar way described by Auer and Mihalas 
(1969) and Mihalas (1978). 


2.2. ATOMIC MODELS 


In the present study, we have assumed that the entire prominence material consists 
solely of hydrogen, sodium atoms and their ions. Hydrogen atoms are the dominant 
contributor of free electrons in prominences, and our main concern is to see in what 
physical conditions the observed emission lines of two atomic species with a large 
difference in ionization potentials and atomic abundances can be explained consistently. 
Therefore, we may be justified in neglecting the existence of other metallic species. 
The hydrogen atomic model used in our work is the four bound levels plus continuum 
representation. For sodium, we assume 7 levels (36, ЗР, 3Р5, 48,5, 3D, АР, 
55 2) and continuum. Transition probabilities of both Н and Na are taken from Wiese 
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et al. (1966), and collisional rate formulae from Sampson and Golden (1970) and other 
necessary data from Mihalas (1967) for hydrogen atom. For the radiative and collisional 
transition rate calculations of sodium, we adopt the same atomic data as listed in 
Landman (1981, 1983). The sodium abundance relative to hydrogen is taken from the 
solar elemental abundances of Lambert and Luck (1978). 

For incident radiations of Læ and Lf, we have used the OSO 8 chromospheric line 
profiles adopted by Vial (1982b) and the total intensity of Ly observed by Vernazza and 
Reeves (1978). The continuum radiation fields have been taken from the tabulation given 
by Heasley et al. (1974), and the incident radiations for sodium lines from the obser- 
vational data given by Pierce and Slaughter (1982). The incident radiation fields are 
assumed to be isotropic and their means are calculated by using a dilution factor, 
ignoring the angular dependence of the radiation. 

Of the numerous allowed transitions, we treat La, LB, Ly, Ha, D1, D2 lines and Lyc 
as optically thick lines and continuum. All the remaining subordinate continua and line 
transitions are assumed to be optically thin. For the optically thin transitions, we regard 
that corresponding incident radiations on the slab surface can penetrate into the 
prominence without any attenuation. Consequently the contributions of those radiations 
to the radiative rate are evaluated as fixed values, with no consideration of transfer 
process. For optically thick lines, we have solved simultaneously the equations of 
radiative transfer, statistical equilibrium and charge-particle conservations to determine 
the distribution of mean radiation field, excitation and ionization equilibria of H and 
Na atoms within the slab of prominence models. 

It is desirable that Lyman lines should be calculated precisely, including the partial 
redistribution and Voigt line profiles. Frequency redistribution effects for La line in 
prominences were pointed out by Milkey et al. (1979) to show that intensity profiles of 
La predicted by complete redistribution and partial redistribution computation differ 
greatly in the wings of the line (see Figure 1 of the paper). But in our calculation, we 
have assumed for simplicity that the lines are formed by a process of complete redistri- 
bution over a Gaussian line profile which is characterized by thermal and microturbulent 
(é = 5 km s^ +) Doppler broadening. 


3. Results 


We have constructed a number of isothermal and isobaric models represented by total 
number density, temperature and thickness as parameters. In our work, 5000 km is 
assumed as a standard thickness of prominence models. The physical properties of the 
models we obtained, are summarized in Table I. N, T, and P denote the total particle 
number per unit volume, temperature and gas pressure of models, respectively. The 
notation (N(i)» and (N,/Ng), „ represent the line-of-sight average of occupation 
numbers in energy state i, and the ratios of proton to neutral hydrogen number Үн at 
the surface and the center. For dense and cool models with a large variation in the 
ionization degree from the surface to the center, the average, (N(i)» is the value in the 
inner part of the slab, since the change in populations is generally considerable around 
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the surface but negligible in the interior. Table II shows the optical properties of 
corresponding models concerning the total optical depth, total intensity of line and 
intensity ratio of D1 to D2. Total emergent intensities of each line are calculated by 
the formal solution of transfer equation: 


оо m 


E- | dv | S,(t,) expC- t,/p) аги, (5) 


where radiations emerging normally from the slab are considered only, by putting the 
projection angle и to 1. 

We find, as shown in Table II, that emissions of all considered lines become stronger, 
in proportion to the rise of the temperature and density. Compared to the hydrogen lines, 
the optical depth of D2 line is always below 1, even in the models with high density. 
Tables also show that the ionization degree of hydrogen at a temperature of 4000 K is 
quite low, resulting in much weaker intensities of NaI D lines at this temperature than 
the smallest observed values of Landman (1981) or Nikaido and Kawaguchi (1983). 
Low pressure models investigated intensively by Heasley and Milkey (1978) yield 
generally low Мат D line intensities, in spite of the high ionization rate. 

We have examined in detail the dense and cool core with total density of 
1 x 10:2 cm °, a temperature of 5000 К and a thickness of 5000 km. Depth depen- 
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Fig. 1. Variation of radiation fields for hydrogen lines and continuum within the slab of a model prominence 
(М = 10'*cm~?, T = 5000 K, D = 5000 km). Abscissa represents the depth of a slab from the surface to 
the center (2500 km) and ordinate, mean intensities denoted by J. 
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Fig. 2. Distribution of populations in the first three energy levels and continuum for the H and Na atoms, 
including the electron number. Model parameters in use are the same as in Figure 1. 


dencies of radiation fields for Lyman lines and continuum including Ha line are shown 
in Figure 1, which indicates that all incident radiations penetrate into a substantial depth 
in the slab. Figure 2 shows the populations of hydrogen and sodium atoms in the first 
three energy levels and continuum. The variation pattern for the first three energy levels 
of sodium atom is the same as that for the free electron populations, implying that 
collisional recombinations of the ionized sodium with electrons are very dominant. As 
seen in the above figures, the most striking features are the very steep gradients in 
radiation fields and populations which occur close to the surface. These steep gradients 
arisefrom the considerable attenuation ofthe Ly continuum radiation to ionize hydrogen 
gas. We notice that sodium atoms are entirely ionized even at a temperature as low as 
5000 K by EUV radiations, while almost all hydrogen atoms remain in the ground state, 
reflecting the large difference in ionization potentials between both atoms. 

Keeping density and thickness of the model unchanged, we made the same analysis 
for a temperature of 8000 K. As shown in Figure 3 and 4, radiation fields and level 
populations change slowly with an inverse gradient compared to the previous low 
temperature model. As mentioned by Heasley and Milkey (1983), this can be understood 
by the fact that the radiation temperature (~ 6500 K) of Ly continuum incident on the 
surface runs close to the internal temperature at the central part of the model where its 
optical depth becomes large. Variations of source functions for Ha and D2 lines are 
plotted in Figure 5 with their contribution functions which are defined by the integrand 


O Kluwer Academic Publishers * Provided by the NASA Astrophysics Data System 


EMITTING REGION OF SODIUM LINES IN SOLAR PROMINENCES 55 


-3.00 
-3.00 


H ALPHA 


-4.80 
-4.80 


О e 

e (fa) 
= 1 LY ALPHA T 
(9 c5 
c) c9 
- © © _J 

=r zT 

eo co 

l t 

© e 

сч а” 

© o 

n т 

о] TOTAL DENSITY = 5 

© | TEMPERATURE (K) = o 

с! | THICKNESS (КМ) = a 

I І 

1x10? 3x10? 1x10! 3x10! 1x10* 3X10 1x10? 


DEPTH (KM) 


Fig. 3. Same as Figure 1 for the model prominence with the parameters of N = 10'? cm ?, T = 8000 K, 
and D = 5000 km. 
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Fig. 4. Same as Figure 2 for the model prominence with the parameters of N = 10'? cm ?, T = 8000 К, 
and D = 5000 km. 


© Kluwer Academic Publishers * Provided by the NASA Astrophysics Data System 


56 KAP-SUNG KIM 


© e 
= сч 
! ! 
TOTAL DENSITY 
TEMPERATURE (К) 
THICKNESS (KM) 
w ou 
Т; СІН ALPHA) e 
C(NR D2) 
5 e zr 
T ~ S (H ALPHA) EIN 
к 1 с 
| сс 
ii D 
а. єз 
"E. о ©? 
X © uc» 
uU» 1 Ы 
al 
C5 S (NA 02) 
© zu 
—! 7. со 
оО Ё 
— со 
[| [| 
(өт) 
. Ce 
(e " 
= T 
D 1000 2000 3000 4000 5000 


DEPTH (КМ) 


Fig. 5. Contribution and source functions for Ha and NaI D2 lines seen at the surface (0 km). Left-hand 
ordinate represents the contribution C defined by the integrand of Equation (5) and right-hand ordinate, 
the source S. 


of Equation (5), S(t)e~ *. As shown in the figure, there are no changes in the source 
functions of Hg and NaI D2 lines. The contribution function for Ha line is large around 
the surface of the model, while that for the optically thin D2 line is practically constant 
throughout the interior. Such a different behavior of the contribution functions for both 
lines is due to the large difference in optical thicknesses of the lines. 

The influence of geometry upon the line emission has been examined by considering 
the fine structure of prominences. That is, we divide the prominence model represented 
by one slab into a number of identical pieces with a given size and apply to each piece 
the same boundary conditions as assumed initially in the model calculations of Tables I 
and II. The radiative interaction between filaments is not taken into account in our work. 
Table III shows the results obtained under the assumption that и filaments with the 
same thickness are lined up in a column along the line-of-sight. D and n denote the 
geometrical thickness of individual filament in kilometers and the number of filaments 
standing in the line of sight direction, respectively. As seen in the table, there is no effect 
of the fine structure upon the ionization rate of hydrogen and the Nar D2 line intensity 
for the model with a temperature of 7000 K, because of the constant distribution of 
populations within the whole interior. However, for the model with a large change in 
the ionization rate from the surface to the center, we can expect an enhanced D2 
intensity by considering the filamentary structure. For instance, the D2 intensity emitted 
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TABLE III 
Effect of geometrical thickness оп NaI D2 line 


EMITTING REGION OF SODIUM LINES IN SOLAR PROMINENCES 


D (km) Т = 5000 К М = 1Е12 Т = 7000 К М = 1Е12 
(Ne? E(D2) (Ne? E(D2) 
5000 1 7.35E9 1.009E2 2.07Е11 2.118E3 
2500 2 7.98Е9 1.096Е2 2.07Е11 2.117E3 
1250 4 1.03E10 1.415E2 2.06E11 2.11453 
625 8 1.50E10 2.060E2 2.06Е11 2.108E3 
313 16 2.45E10 3.363E2 2.05Е11 2.099E3 
156 32 3.01Е10 4.130Е2 2.03E11 2.083E3 


E(D2) is a total intensity of NaI D2 line emitted from n filaments with a thickness D standing in the 
line-of-sight direction. 

< №. > colums shows the change in the line-of-sight average of free electrons with the geometrical thickness 
of model filaments. 


from a group of filaments (T = 5000 K, D = 156 km, and n = 32) is about 4 times higher 
than that from one slab model of D = 5000 km, because the number of free electrons 
from hydrogen atoms in a filament is highly increased by the geometrically thin structure. 

We can make the Ly continuum radiation penetrate the filament more freely, by 
reducing its density. Figure 6 shows the distribution of radiation fields for a filament with 
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Fig.6. Same as Figurel for the model prominence with the parameters of N = 3 x 10!! cm ^ ?, 
T = 4000 K, and D = 300 km. 
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a thickness of 300 km, a temperature of 4000 K and a density of 3 x 10!! cm ^? lower 
than in previous models. It is remarkable that the incident radiations pass through the 
whole interior without a large attenuation, for this model with a temperature as low as 
4000 К. Although we can raise the ionization rate by reducing its density, strong NaI 
D line intensities are difficult to obtain. This is caused by the fact that the number of 
sodium atoms is reduced in proportion to the decrease of the total density. 

In addition to the above mentioned facts, we find that the source function for Ha line 
is approximately constant, regardless of the density and the temperature. This property 
of constant source function for Ha is consistent with the result of Fontenla and Rovira 
(1985), and can be used very conveniently for the empirical analysis of the observed 
spectral line. For NaI D2 line, its optical thickness is so thin that the effect of radiative 
transfer may be negligible, except for the model case of extremely high density. 
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Fig. 7. Plot of observed Nai D2 integrated intensity against that for Hf line. Full lines are the theoretical 

sequences calculated by models with a constant density of 3 x 10!!, 10'2, and 3 x 10? cm ^?, and variable 

temperatures. Symbol ж marked on the curved lines by 1, 2, and 3 represents the temperature of 8000, 7000, 
and 6000 K, respectively. 


We have compared the results of our calculation with the observations. A plot of NaI 
D2 integrated intensity against that for Hf line is shown in Figure 7. Symbols of О, A, 
+, X, and <> in the figure represent the spectral data observed by Sergeeva and 
Yakovkin (1964), Sergeeva (1968), Yakovkin and Zeldina (1974), Yeh (1961), and 
Kubota and Kureizumi (1986). respectively. We detect the observed proportionality 
between the NaI D2 line strength and the intensity of Hf, as already found in Table II. 
We have attempted to match the observations with some of our models. Full lines in 
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Figure 7 are the theoretical sequences calculated by models with a constant density of 
3 x 1011, 10!? and 3 x 10!? ст 3, and variable temperatures. Calculations for models 
with a density higher than 10!? cm ^? are in good agreement with the observations. On 
the contrary, low density models cannot fit well the observations, unless their thick- 
nesses are assumed to be unreasonably large. 

The electron density obtained from our models with a temperature higher than 6000 K 
generally lies in the range 10!9—4 x 10!! cm-?, except for the model with 
N = 3 x 10? cm ^? and T = 8000 К. This is nearly consistent with the result of Chultem 
and Yakovkin (1974) and Hirayama (1985). Much lower electron densities of 
N, ғ 7 х 10? cm 3 have been derived by Sahal-Brechot (1984) from the Hanle effect, 
but our models with a low electron density do not reproduce well the emissions of NaI 
D lines. Our ionization degree varies in the wide range 0.9-2 x 10-2. Landman’s value 
for the rate (0.16) corresponds roughly to the model with N= 10!! cm ^? and 
T = 6000 К or N = 3 x 101? cm ^? and T = 7000 К. For models with density higher 
than 1012 cm ^? and temperature above 7000 K, the calculated optical depth of Hf 
exceeds 1 and t(Ha) becomes very large, compared with observed values (10-60, 
Hirayama, 1964; 7-32, Yakovkin and Zeldina, 1974). Our calculated values of the 
intensity ratio of D1 to D2 line (~ 0.62) agree very well with the observed ratios shown 
in the work of Nikaido and Kawaguchi (1983). 


4. Discussion 


In our model calculations, we have assumed 7 bound levels plus continuum for the 
sodium atom. The truncation effects of remaining upper levels were pointed out by 
Landman (1986). He showed that the neglect of the upper level group (> 5P) in Na atom 
leads to roughly a factor of 2 underestimate in the evaluation of lower level populations. 
We admit that our calculated results on the total intensities of Мат D lines may be 
underestimated by a factor of about 2. 

We will discuss the role of some physical parameters related closely to the sodium 
lines formation. The strength of sodium lines should be proportional primarily to the 
total amounts of sodium atoms along the line-of-sight in the excited states of 3P,/ 3/2 
levels. Consequently, the total number density and the thickness of the slab are the major 
factors to control the intensities of Мат D lines. 

In the statistical equilibrium solution, atomic populations are determined by the 
temperature, radiation fields and electron density. For a fixed number of total particles 
and electrons, the strength of Nar D lines increases as the temperature decreases, due 
to the corresponding rise of the populations in ЗР, 3/2, including 35, level. In view 
of this fact, it is natural that Nikaido and Kawaguchi (1983) derive an emitting region 
with a low temperature (~ 4000 K) for strong Nar D lines, on the ground that they 
constructed the prominence model with a fixed electron density (N, = 0.1N). Actually, 
however, free electrons vanish, in proportion to the decrease of the temperature, as 
shown in Table II. At the same time, emission of sodium lines falls down, as the number 
of electrons decreases with temperature. 
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We would like to emphasize here the role of free electrons in the formation of sodium 
D lines. Electrons, which are mainly supplied from hydrogen atoms ionized by Lyc 
radiation, are again recombined with ionized sodium, throughout dominant collisional 
recombination processes. As pointed out by Landman (1983), the magnitude of electron 
density yields the large populations in ЗР; з» levels including the ground level, 
according to the relation of V(3P;) ~ Ne N(Na* ) depending slightly on the temperature. 
Thus if free electrons (or ionized sodium) are numerous, we can expect strong emissions 
in NaI D lines. Then, the number of electrons should be determined by the simultaneous 
equations of (1), (2), (3), and (4) in our work. Actual electron density calculated in this 
way has a strong dependence on the model parameters such as temperature and density, 
as shown in Table I. It is apparent that the ionization degree adopted in previous work 
(Nikaido and Kawaguchi, 1983) as a parametric independent value, should deviate 
largely from the real one for models of temperature below 7000 K. 

This can be proved more evidently by the inspection of the variations of the ionization 
rate with temperature. Figure 8 shows the behavior of the ionization rate obtained by 
the non-LTE calculation of hydrogen. The rate is almost constant with a value of about 
0.3 for the model with 7000 K, but at 4000 K the rate is so small that the total number 
of electrons supplied from hydrogen, even including the contribution from sodium are 
no more than 10 ^? times of the total number of hydrogen atoms. Actually, great care 
must be paid when one estimates initially the electron density as a fixed value in 
non-LTE calculation of heavy elements. 
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Fig.8. Ionizationratio, N(H * )/N(H1)ofamodel prominence(N = 10!? cm~ ?,D = 5000 km)forthe various 
temperature values. N(H + Jand N(H 1)denotethe number ofionized hydrogen(protons) and neutral hydrogen, 
respectively. 
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Provided that free electrons are abundant, it is expected that Nar D lines could be 
formed sufficiently, even in prominences with a temperature below 5000 K. We have 
attempted to increase the electron density by reducing the thickness of the slab. The 
expected increase of the ionization degree was obtained for a cool and dense model with 
a large change in radiation fields and populations from surface to center. Aside from 
changing the thickness, the decrease of total density also leads to the fair rise of the rate, 
due to the easier penetration of Lyc radiation. But the variation in strength of NaI 
D lines is insignificant, since total amounts of sodium atoms are depleted by lowering 
the total number density. 

The strength of the EUV continuum radiation has a great influence on the ionization 
equilibrium of atoms. In our work, we used the continuum radiation field observed in 
quiet Sun. Solar EUV(A < 912 A) and UV(A < 2414 A) observational data show that 
the change of their continua from quiet regions to coronal holes is very small. However, 
we find that the Ly continuum intensity above active regions comes up to about 8 times 
of that above quiet ones, as shown in Vernazza and Reeves (1978), while the UV 
continuum (A < 2414 A) varies slightly throughout the two different regions. According 
to the calculation with the changed Ly continuum data, the number of free electrons 
from hydrogen is increased by a factor under 2 at the surface. For the most part of the 
interior, there is no change in the number of free electrons, since the variation of Ly 
continuum occurs only around the surface. Consequently, even if the continuum 
radiation data from an active region are used, the change in the line-of-sight average of 
free electron number < №, У and the emission strength of Nar D lines is insignificant, 
compared to the quiet region. 

Emergent intensities of Lyman lines wholly originate from the narrow sheath of 
prominences, due to the extremely large optical depths of those lines. In fact, our uniform 
models with a temperature lower than 8000 K cannot produce the central reversal of 
La and Lf lines confirmed by satellite observations (Vial, 1982a), since the temperature 
in the vicinity of the surface is uniformly fixed to that of the cool core. If a realistic model 
structure with a steep temperature rise near the surface is considered, one can reproduce 
well the centrally reversed profile. Even for uniform models, reversed profiles can be 
obtained, if temperatures higher than 8000 K are assumed. Anyway, we would like to 
stress that our main concern is not a hot sheath where the Ly lines are formed but a 
cool core of the prominence to examine the emission of sodium lines. 

In any case, numbers of prominence models with strong emission of Мат D lines are 
possible by the suitable combination of physical parameters. Among the numerous 
models supposed, we try to determine reasonable physical conditions for the NaI 
D lines emitting region, by using the observational constraints of Ha and NaI D lines 
in quiescent prominences. In fact, far more models must be simulated to examine the 
physical conditions which can reproduce well the observed quantities. Even though a 
sufficient number of examples has not been presented in our work, we can make the 
contour map which represents the iso-intensity of sodium D2 line and the iso-optical 
thickness of Но line, by an interpolation based on our results of model calculations. 

Figure 9 shows the diagnostic diagram on the physical conditions of the emitting 
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and dashed lines represent the iso-contours of optical depths for Ha, intensities for Мат D2 line and gas 
pressures, respectively. 


1X10 


region in model prominences with the fixed thickness of 5000 km. Thin lines in the figure 
represent the contours for iso-optical depth of Ha with value of 1 and 50. The band 
formed by those lines stand for the physical conditions of Ha line emitting region in 
prominences, which belong to the range of the observational constraints for the line. 
Contours of iso-intensity of Мат D2 line with 6 x 107, and 10^ ergs cm~? 57! ster ^! 
are drawn in the figure with thick lines. Since the intensity of Мат D2 line observed in 
quiescent prominences lies in the range 6 x 10? — 10* ergs сп” 2657! ster ^ !, the 
emitting region of sodium lines in quiescent prominences should be located in the zone 
enclosed by the contour lines. We also carry out the check of the gas pressure of 
prominence models. As drawn with dashed lines in the figure, the isobaric line corre- 
sponding to the coronal pressure value of 0.18 dyn cm ^? lies outside the Мат D2 line 
emitting region. It is certain that Мат D lines with an appreciable strength should be 
emitted in the high pressure region on the diagnostic diagram. Referring to the results 
of Landman (1983, 1984), we restrict the upper boundary of the pressure for our models 
to 4dyn cm "2, As shown in Figure 9, physical conditions of Но and Мат D lines 
emitting regions, which satisfy both optical and physical constraints in quiescent 
prominences, are situated within the range represented by the shaded area. Small filled 
area represents the physical conditions where quiescent prominences can emit the strong 
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Nal D lines with an intensity over 10* ergs cm? s7! ster ^ !. In fact, such a strong 


emission of D lines may be infrequent in quiescent prominences. 

From the above diagnostic diagram, we can estimate that sodium D lines are difficult 
to be formed in the cool core of prominence with a temperature below 5000 K. In order 
that the cool part of the prominence with a pressure below 4 dyn cm ^? emit the sodium 
D lines with an intensity between 6 x 10? and 10* ergs cm 2687! ster ^ |, along with a 
Ha optical depth between 1 and 50, a density higher than 4 x 10!! cm ^? and а 
temperature over 5000 K are necessary at least. Our conclusions on the emitting region 
of sodium lines are supported implicitly by Landman's results (1983, 1984) for the 
principal plasma parameters derived from the empirical analysis of optically thin lines 
including Мат D lines. 
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